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Abstract 

A renormalization-group and bosonization approach for a multi-band Hub- 
bard Hamiltonian in one dimension is described. Based on the hmit of many 
bands, it is argued that this Hamiltonian with bare repulsive electron-electron 
interactions is scaled under specific conditions to a model in which supercon- 
ducting fluctuations dominate. 
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The Hubbard model plays a paradigmatic role in understanding many-body correlation 
effects. Stimulated by the discovery of high-Tc superconductivity, it has been subject of 
renewed interest and its multi-band versions with purely repulsive interactions have been 
examined as possible realizations of superconductivity [|l|]. However, the first study of multi- 
band Hubbard models date back to Van Vleck who evoked intra-atomic interactions 
between orbitally degenerate states as an explanation for ferromagnetism. In general, these 
models are fascinating because, due to orbital degeneracy, they incorporate a variety of 
intriguing physical phenomena, including, for example, the coexistence of orbital superlattice 
and magnetic orderings 0. One potential application of multi-band models could be doped 
fuUerenes — e.g.. Ceo with three- and five- fold degenerate tiu and hu bands — where the 
consequences of band degeneracy are certainly important in understanding the physics of 
these novel materials. 

While the above studies mostly consider narrow bands or only a few degenerate bands 
within mean-field approximation, Muttalib and Emery |Q studied recently a two-band model 
of spinless electrons in one dimension by abelian bosonization and renormalization-group 
analysis. Interestingly, they found that superconductivity can exist with repulsive electron- 
electron interactions. They also considered the question of decoupling of dynamical degrees 
of freedom which they found to occur only at some specific values of coupling constants. 

Motivated by the above examples, we focus in this paper on a novel, degenerate multi- 
band Hubbard model for fermions with spin degrees of freedom. We first derive one-loop 
renormalization-group equations and analyze their behaviour. Solvable limits of this model 
are then studied by bosonization. We argue that, for a particular choice of parameters, 
the model exhibits enhanced superconducting fluctuations. We conclude by noting that the 
decoupling of dynamical degrees of freedom does not occur in general, in contrast to the 
one-band Hubbard model. 

Specifically, our starting point is the Hamiltonian 
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where Cnaa is the fermion operator for an electron of spin a with band index a and pna = 
So- c]^ao■<^nao■ at site 71. We associate the band degeneracy to colour degrees of freedom 
and spin degeneracy to flavour degrees of freedom; the corresponding indices take values 
a = 1, . . . , Nc and a — 1, . . . , Np. Ua/s is the on-site electron-electron interaction which 
takes into account possible couplings between various bands, and t is the hopping amplitude 
which has been assumed to be the same for all the bands. In the absence of electron-electron 
interactions, the model has the U(7Vc x Np) symmetry which is preserved for Uap = U. 
By including more than one band and allowing the broken colour symmetry, interesting 
possibilities arise, one of which is enchanced superconducting fluctuations, as seen below. 
In the following, we always assume the full spin-rotational symmetry. 

Our primary results can be qualitatively explained by considering the Hubbard interac- 
tion in the U(A^f) -invariant form of 

W± E ( EP-)' + - EP-Pna, (2) 

n a na 

where we have assumed that the diagonal interactions can be parametrized by U\\ and the 
off-diagonal interactions by U±. Subsequently, we will be mostly interested in cases where 
U±> while may have an arbitrary sign. Our results are: (i) In the limit of Nc —>■ oo, 
U± scales to zero while U\\—U± is invariant under the scaling. Thus, the model will scale to 
Nc independent systems whose on-site electron-electron interactions are either attractive or 
repulsive, depending on the sign of U\\ — U±. This is expected, because the mean-field ap- 
proach should work in this limit and so the first term in Eq. (2) is qualitatively unimportant. 
(ii) For finite but large Nc, the system again scales so that U± decreases and renormalized 
diagonal interaction may be positive or negative. However, if initially U\\ < U±, there exists 
a limiting value for the effective U± after which further scahng causes it to increase. This 
critical value scales as l/Nc- 

Current Algebra and the Action. We will use the current-algebra formalism in the calcu- 
lations because it provides a compact notation and allows operator-product expansions to 
be used efficiently in deriving the scahng equations as short-distance degrees of freedom are 
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integrated out. Since we are interested in low-energy and long- wavelength phenomena, the 
spectrum is linearized at the Fermi energy and the fermionic degrees of freedom are expressed 
by slowly varying fields ipaa{x±), where the left (+) and right (— ) moving electrons are la- 
beled according to their arguments, x± = t ± x. We define Jaa''ix±) = ■'4^tai^±)i^a' a' ix±): 
which present the left and right moving currents; the upper (lower) indices of the current 
refer to the flavour (colour) degrees of freedom. The colons denote normal orderings with 
respect to the filled Fermi sea of the noninteracting system. Flavour and colour currents are 
then defined naturally as J"(x±) = Taa"^aa i^±) J'^i^i) = T^a"^aa'i^±)^ respectively 
(hereafter a sum over repeated indices is implied) . Hermitian matrices and form a Lie 
algebra of SU(A''i?) and SU(Ai'(7) groups. Note that it is convenient to add the identity ma- 
trix to the above sets of generators; for example r° = 1/ y/Np . Therefore, both currents 
actually give the U(l) current with a = A = 0, whereas J" and are the true SU(A'^f) 
and SU(Ai'c) currents for nonzero a and A. It is easy to compute the current commutation 
relations, equivalent to the Kac-Moody algebra. However, more useful is the corresponding 
operator-product expansion, 

271 {z± — w±) 

This follows from the commutator and can be shown by means of the conformal symmetry. 
Here we assume radial quantization where the spatial dimension is compactified by mapping 
the space and time coordinates to a cylinder of radius 27r; z± = exp(t ± ix). Below, we 
continue to use the fiat space-time coordinates for the action but do the actual calculations 
in the complex plane. 

The free part of the action is 

= .r Iat E / d'x[r{x,)r{x,) + J^{x,)J^{x,)], (4) 
and the action describing the electron-electron interaction is 



4 



Si — J (fx (^—U''^pJ'^p{x+)Jp'^{xJ) + uj^f^J^^{x+)Jp^{x-) 

+ lf45[:^L(^+)V^«.(^-): + h.c] (5) 

where (fx = dxdt and the units are defined so that the Fermi velocity becomes hvp = 1. In 
contrast to the Hubbard model, we will allow the coupling constants U'^^^ to be independent, 
as implied by various physically important interactions — for example, by the nearest- 
neighbour electron-electron interaction. 

Renormalization- Group Equations and the Flow. We follow the usual procedure in which 
the partition function, Z = J V[ilj]e~^ , is regularized by introducing a ultraviolet cutoff i for 
short distances. As the short-distance degrees of freedom are integrated out continuously 
by using the operator-product expansion (|^), the coupling constants are correspondingly 
renormalized so that Z remains unchanged This procedure then gives the one- loop 

scaling equations for the couphng constants: 

dm =-i9!± + 9L) + [{9f±+9l±)-9n]/Nc 

digii. = -{gj± + 9i±) + [{g2\\ - 2^i|| - g2± + 2gi±)gi± + {2gs± - g3\\)g3±]/Nc 

dm = {gl^^-g%)/2Nc 

dm = {gl^-gl^)/2Nc 

dm\ = -25'i±fl'3± + [{'^g2\\- gi\\)gn + '^giLgzi]/Nc 

dig3± = -2gi±g3± + [{g2\\ + g2± - 2gin + 45(i_l)5'3± - fl'3||5'i±]/^c, 
where gk\\,± = NcU'^\]_/'k (/c = 1, . . . ,4) and Np = 2 is chosen. The coupling constants are 
functions of the ultraviolet cutoff 9^ = d/d\ogl. As usual, gia is the backward-scattering 
constant, g2a is the forward-scattering constant, and g^a is the Umklapp-scattering constant 
(a = ||,_L). Because g^a corresponds to the scattering processes which violate momentum 
conservation by 4fc^, it is important only if Akp is equal to the reciprocal lattice constant; kp 
is the Fermi wavevector which has been assumed to be the same for all the bands. We have 
ignored g^^a because its only effect is to renormalize Fermi velocities; although this could be 
taken into account, 5'4q does not enter at one-loop order since / dz/ = 0. 



Because the parameter space of the system is very large, we have investigated the scahng 
behaviour of the system in two special cases: 

(1) First, we consider the system with the broken orbital U{Nc) symmetry so that the 
action is only U{Nf) invariant. We further set the forward-scattering coupling temporally to 
zero, g2\\,± = 0, because the backward- and Umklapp-scattering coupling constants are the 
most relevant ones in terms of that their scaling behaviour is nontrivial. By letting initially 
g\\,± = gk\\,i. {k = 1, 3), and 5'2||,± = 0, these properties are conserved by the renormalization- 
group flow. Thus, 

degii = -2gl + i2gl-g^^)/Nc, 

dig± = -2gl + (Ag^ - ?,g\{)gjNc] 
with Np = 2. It is easy to see from these equations that the renormalization-group flow 
may change the sign of g\\ but not the sign of g^^. In fact, g± = ^ forms a fixed point 
for this coupling constant. This, and the line g\_ = g\\ act as separatices: (i) First, let 
fi'x > fi'll > 0. Now there is also another relevant solution for the equation dig± = 0: 
g± = —3g\\/{2Nc — 4). This solution yields a turning point for the scaling trajectories at 
which the initially decreasing gj_ starts to increase. For large Nq, the turning value of g± 
scales as I/Nq- Note that g\\ decreases monotonically and, at the turning point, gi\ < 0. (ii) 
In the region g\\ > g± > 0, the system scales towards a noninteracting model {g\\,± = 0), 
corresponding to the Nq free fermion theories, {iii) Finally, g± < leads always to the 
strongly attractive interactions: g\\ and g± decrease without bound at one-loop order. 

To illustrate the general behaviour of Eqs. (|^, we integrate them numerically and plot 
the renormalization-group flow for Nq = 10; see Fig. 1. For g± > g\\ > 0, it clearly shows the 
tendency of the renormalization-group flow to lead to an effective Hamiltonian which con- 
tains weakly coupled one-band Hamiltonians with intraband electron-electron attractions. 
It is then convenient to diagonalize each of these Hamiltonians separately and to take into 
account the off-diagonal interactions perturbatively. However, this assumption relies on the 
fact that the "unperturbed" system has a gap the perturbation expansion to work. These 
results apply to Eqs. (||) whose qualitative behaviour is similar for finite g2\\,±- 
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(2) We next consider the model with the full \J{Nc x 2) symmetry: t/^^ = U^''\ k = 
1, . . . ,4. In this case, the interaction term becomes 

Sj = J (px[ - u'-^^r^{x+)r^{x^) + f/(2) j(x+) 

+ [J{x+)J{x+) + J(x_) J(x_)] (8) 
+ Umklapp^, 

where J"^(x±) = 'r^a''^aa'Jaa'i^±) is the colour- flavour current and J = is the U(l) 
current. As an example, we assume that the bare interactions are repulsive {gk > 0). In 
this case and at half filling, gi can scale to large negative values, g2 is renormalized only 
weakly, and g^ decreases and starts to increase approximately when gi becomes negative. 
Note that the renormalized value of g^ is always positive. It is very difficult to draw any 
conclusions, even at gi = for g^ (the half- filled case), although the interesting feature is 
that the backward-scattering coupling constant gi scales towards large negative values which 
may suggest enhanced density-wave fluctuations in colour or flavour channel. In general, 
an analytical solution does not seem to be possible when gi and g^ are the most dominant 
coupling constants. 

Bosonization. We now turn to a limit where the backward and Umklapp terms are ei- 
ther zero or scale to zero, so that the multi-band model is exactly solvable by bosonization 
I^J^. In particular, we consider the U(A'^p) -invariant model where the '[J{Nc) symmetry is 
broken by letting U\\ ^ U±. For our purposes, abelian bosonization is sufficient. At zero 
temperature, the most dominant instabilities are either 2kF density-wave (DW) or super- 
conducting (SC) instabilities whose correlation functions have power-law singularities 
at low energies and small momenta with exponents 

/iDw = 2(1 - 69) /Nf + 2(69- 9)/NcNf, (9a) 
/isc = 2(1 - 1/69)/Nf + 2(1/69 - 1/9)/NcNf, (9b) 

where 69'^ = (1 + 6U4 — 6u2)/{l + 6U4 + 6U2) and 9"^ = [1 + U4 — U2)/{1 + U4 + U2)- We 
have defined 6uk = Nf{uI''^ - uP)/A7r and Uk = Nf[uI''^ + {Nc - l)uP]/47r, {k = 2,4). 
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In these expressions, ?7||^| are the renormahzed values of the forward-scattering constants. 
Clearly density-wave fluctuations always dominate over superconducting fluctuations for 
f/J^') = U^l'^ > and vice versa. However, if < U^^ < U^l\ superconducting fluctuations 
win, provided that Nc > {66/6 — I)/ {66 — 1) and U^^ < where the latter condition is 
needed to guarantee a sensible large- Ai'c limit. 

In contrast, if the large- A^c limit is taken so that gk\\.± = NcNfUI^^]_/ 2tx are kept constant, 
the condition, < g2\\ < g2L, is no longer sufficient to have diverging superconducting 
correlation functions at low energies and small momenta. As an example, let g\\^± = gk\\,i. 
{k = 2, 4) and consider the weak-coupling limit: g\\^± ^ 1. We find that the superconducting 
correlation functions are the most singular ones for the repulsive interactions and Nc ^ 1, 
if the condition g\\ < gj_/2 + 0{gj_) is satisfied |]rT |. 

It is interesting to note the connection between the above consideration and the 
renormalization-group approach in the case of the repulsive interactions: they both im- 
ply that the necessary condition for superconducting fluctuations to dominate is f/y < U± 
and Nc to be large enough. 

Colour-Flavour Separation. The scaling equations show no apparent separation into in- 
dependent sectors. This not so surprising if we consider the model with the U{Nc x Np) 
symmetry. This coupling of the colour and flavour degrees of freedom — a manifestation 
of the fact that the colour and flavour currents are not independent dynamical variables — 
follows because the colour-flavour current commutes neither with the colour nor with the 
flavour current [1^. However, away from half fllling, the U(l) current describes an inde- 
pendent dynamical degree of freedom; J commutes with the other currents. The coefficient 
multipling the J J term in Eq. (^) is U^^^ = U^'^^ — U^^^/NcNp. Consequently, in the one-loop 
renormalization-group equations (^, (yfc = NcNpU'^^'> 1271 is indeed decoupled from the other 
degrees of freedom. Moreover, it is a scaling invariant. In the one-band Hubbard model, the 
spin-charge separation is evident even at half-fllling, because the Umklapp term is a chiral 



^\]{Nf) singlet |T3 
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Relation to the Muttalih-Emery Approach Our approach differs in two respects 
from that of Ref. [Q: we consider an arbitrary number of bands with the spin degrees 
of freedom but ignore the interband Umklapp scattering. The interband Umklapp scat- 
tering would allow versatile possibilities; for example, one could include terms such as 
U^p Jd'^x J^p{x+) J^p{x J) . Since the scaling equations have to be invariant under the phase 
transformation ifjaa — > e^^°"^acTi = is a fixed point and the other scaling equations 
contain only even powers of U^^p. The question remains whether these coupling constants 
are (marginally) relevant. Because they do not originally appear in the Hubbard model, we 
have assumed that their bare values are negligible and that they do not change the scaling 
behaviour qualitatively. On the other hand, in the model considered in Ref. [Q, they were 
found to be important. Our motivation was to illustrate that even the Hubbard model can 
show enhanced superconducting fluctuations if the spin and multiple bands are assumed. 

Conclusion. We have shown that the existence of multiple bands can lead to enhanced 
superconducting fluctuations. Based on the one-loop renormalization-group analysis, we 
found that the renormalized diagonal electron-electron interactions are attractive and the 
off-diagonal electron-electron interactions are of order 0{N^^), if bare values of the coupling 
constants are such that U\\ < U±. This is consistent with the result obtained by studying 
the Hamiltonian in the exactly solvable limit. 

A.V.B. acknowledges the support of J.R. Oppenheimer fellowship. This work was sup- 
ported by the U.S. Department of Energy. 
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FIGURES 



FIG.l Scaling trajectories in the {g\\, g±) -plane for the half-filled, Nc-hand Hamiltonian for 
Nc = 10, as given by Eq. (0) {Np = 2). The equation g± = —3g\\/{2Nc — 4) is shown by a 
dashed line. 
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